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Abstract 

We study the equilibrium shapes of vesicles, with an in-plane nematic order, using a Monte-Carlo scheme and 
show that highly curved shapes, like tubes and discs, with a striking similarity to the structures engendered by 
certain curvature sensing peripheral membrane proteins, can be spontaneously generated by anisotropic directional 
curvature with nematic disclinations playing and important role. We show that the coupling between nematic 
order and local curvature could lead to like defects moving towards each other and unlike defects moving away, in 
turn leading to tube formation. Thermally induced defect pair production lead to branched tubular structures. 
It is also shown that helical arrangement of the membrane tubes, with nematic field spiraling around it, is a 
dominant soft mode of the system. 

PACS numbers: PACS-87.16.D-, Membranes, bilayers and vesicles. PACS-05.40.-a Fluctuation phenomena, random pro- 
cesses, noise and Brownian motion. PACS-05.70.Np Interfaces and surface thermodynamics 



* Electronic address: ram@physics.iitm.ac.in 
^Electronic address: sunil@physics.iitm.ac.in 
1-Electronic address: ipsen@memphys.sdu.dk 



Fluid interfaces with an inplane orientational order can sustain amazingly complex morphologies. Apart 
from the interesting physics that it can offer, in terms of disclination dynamics on deformable surfaces, 
study of such interfaces are also contribute to understanding morphologies of cellular organelles. 

The organelles of a biological cell have membranes with highly curved edges and tubes, as seen in 
the Endoplasmic reticulum, the Golgi and the inner membrane of mitochondria. Tribulation has also 
been observed, in vitro, in self assembled systems of pure lipids pQ. It has been shown that, macro- 
molecules, which constitutes and decorates the membrane surface, strongly influence the morphology of 
membranes. For instance proteins from the dynamin superfamily are known to pull out membrane tubes 
while oligomerizing themselves into a helical coat along the tube [2J. The BAR domain containing proteins 
in general can induce a wide spectrum of membrane shapes ranging from protrusions to invaginations 
depending on the geometry and interaction strength of the BAR domain [3H5] . 

The standard Helfrich model[6] for membranes, based on mean curvature energy, cannot explain the 
stability of such highly curved structures. Existence of anisotropic bending energy will be the minimal 
requirement to explain the stability of tubular shapes, which could arise from an in-plane orientational 
field on the membrane [7J [8] . Planar orientational order could be intrinsic to the membrane [9] , due 
to the structural properties of its constituents, or could arise as a result of membrane interactions with 
external agents [3_H5]. Lipid molecules tilted with respect to the layer normal, as seen in the L [TO] 

and Pp phases, is a well known case of intrinsic in-plane orientational order. Anisotropic peripheral 
proteins (curvactants), inducing membrane deformations |I2j either by scaffolding or by insertion of their 
amphipathic helices into the membrane [3J, are examples for extrinsic sources of orientational order. 
In general the order can be represented by a p-atic surface field, invariant under rotations by 2ir/p. 
Equilibrium shapes of a surface with p-atic in plane order, have been investigated earlier [131 114) . Mean 
field phenomenological models of orientation order on membrane surface, in the context of lipid tilt and 
chirality, have predicted the stability of tubular membranes and helical ribbons [T5HT5] , More recently, a 
mesoscopic particle model was used to show that presence of BAR domains on membrane surfaces can 
lead to tubulation and vesiculation HH [20] . 

The coupling of curvature to orientational order leads to interesting defect dynamics and thermally 
excited shapes. Thus in addition to its biological relevance, membranes with in plane order also provides 
a unique medium to study disclination dynamics on deformable surfaces. So far such studies have been 
carried out only on surfaces with prescribed curvature |21j . Apart from the interesting physics that it 
can generate these investigations may also provide new routes to new functional materials |22j. 

In this paper we use the Monte Carlo model developed by us [23J [H] to explore the spectrum of 
vesicle shapes engendered by an in-plane nematic field m. The technique developed by us enables one 
to go beyond the mean field regime and also to explore the full nonlinear regime. We scale all energies 
by fcgT. A moderate mean curvature stiffness k = 10 is chosen to enable study of thermally induced 
shape fluctuations. Higher values of k are chosen to study equilibrium shapes in the absence of thermal 
fluctuations. We choose the directional curvature moduli, along and perpendicular to the orientational 
field, as K\\ = 5 and k± — respectively. The results with k± ^ is qualitatively same and is not discussed 

separately. We vary the spontaneous curvature along the in-plane field Cq, to explore different shapes. 
The in-plane nematic field, distributed uniformly over the entire surface, is constrained to be in the ordered 
phase by setting the Lebwohl Lasher coupling constant cll = 3fcgT 25j. For this values of the coupling 
constant, with Ku = 0, the nematic field orients itself on a tetrahedron with four +1/2 disclinations, 
as predicted in [T3]. We will first investigate the effect of the directional spontaneous curvature on 
the equilibrium membrane conformations when the k value is high so that thermal undulations can be 
neglected. 

At high k thermal undulation does not change the shapes significantly and the shapes obtained from 
the equilibrium simulations will be the same as that from a mean field calculation. It is known that 
dynamic triangulation Monte Carlo (DTMC) models, based on Helfrich Hamiltonian 6J, do reproduce 
the mean field phase diagram of vesicles [26 . Here we determine the equilibrium shapes of decorated 
vesicles with an in-plane orientational order, also taking thermal fluctuations into account . As can be 
seen from Figjl] we see tubular and disc like shapes emerging, even when there are no constraints on 
the area to volume ratio of the vesicle. Directional spontaneous curvature and interaction between the 
disclinations are two important factors that determine the shape of these decorated vesicles. Since the 
membrane is self avoiding the constraint imposed by the fixed topology also plays a role. 

Four curvatures can be identified, the curvatures along the principal directions c\,C2 with c\ > ci and 
the directional curvatures, c'^c 1 , along and perpendicular to m. Figl2| shows the distribution P(cj), of 
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(a) Prolate (b) Tubular (c) Disc (d) Caveola like 




(e) Tubular (f) Spirals (g) Narrow neck (h) Broad neck 



FIG. 1: Membrane conformations for different values of C with k = 10, Kn = 5, kj_ = and en, = 3. ( 
See supplementary materials for high resolution figures of nematic field arrangement.) 
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FIG. 2: (Color online) (a,b) Distribution of the maximum and minimum principal curvatures, c\ and C2 

respectively, for different values of c\ . 



Ci=i,2: on a vesicle, for different C values when k± = 0. It is clear from FigMa & b) that for C n > the 

distributions have a single peak, with -P(ci), peaked close to, but less than, the value of Cq, while P(C2) 
is peaked at zero. The resultant value of the principal curvature arising from the competition between 
the resistance to bending imposed by k and the directional spontaneous curvature set by Cq , is shown in 
FigJ^J When k decreases the peak position of P(c\) moves toward that of Cq and then to a value higher 

than Cq for k ~ 0. Note that for k = there are no competing elastic forces and the desired directional 
curvature can be achieved by setting m at an angle <p ^ with respect to the direction of c±. For a tube, 
a non zero value of tp implies a configuration with nematic spiraling around the tube. The angle ip is 
shown as inset in Fig. [3] As expected this angle is a decreasing function of the bending modulus k and 
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saturates to zero for large bending rigidity. 

Negative spontaneous directional curvature (Cq < 0) induces tubes to form into the vesicle. However, 
restrictions imposed by the membrane self avoidance and topology prevent a complete inversion of the 
membrane through inward tubulation. This leads to a bimodal distribution of c 2 , with a peak at c 2 = 
in addition to the expected peak at c 2 ~ Cq (see Figj2] a) . This bimodal peak arises due to presence of 
inward tubes (with most of the vertices with c 2 < 0, c\ = ) and an outward curved surface (with most 
vertices having c\ > 0, c 2 = 0). 




FIG. 3: (Color online) Comparison of the C\ (filled symbols) and directional curvature c" (open 
symbols) distribution, at c\ — 0.4, kj_ — 0, for different values of K. Inset shows the change in the angle 
ip as a function of k. Dotted vertical line marks ci = C = 0.4. 



Another important factor that determines the shape of a vesicle is the positioning of disclinations. We 
know that, for a nematic field on a closed surface with spherical topology, the total disclination strength 
should be 2. There is no such topological restriction on the number of defects themselves and thermal 
fluctuation could excite oppositely charged defect pairs. 
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FIG. 4: (Color online) Distribution of the nematic defect separation (£) for n — 10, £ll = 3 and k± = 0. 



In the absence of any stiffness-orientation coupling, the energy of a pair of defects, on a surface of 
spherical topology, decreases logarithmically with separation [27| and is repulsive when their charges are 
of the same sign [13]. We therefore expect, at low temperatures, to see four +1/2 disclinations at equal 
distance from each other. In Figj4] this is evidenced in the single peak in the measured distribution 
of geodesic distances (-P(O) between the defects. Directional bending stiffness ( n\\ ^ ) alters this 
distribution considerably, since curvature influences the interaction between defects [23 [29] . We see that 
two of the +1/2 defects move towards each other to form pairs as the vesicle deforms into a tube, resulting 
in a peak in the £ distribution at short separation as shown in Fig[4[ This peak moves towards the left 
I 

as Cq increases. The pairs, which now has a total strength of +1, moves away from each other, resulting 
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FIG. 5: (Color online) Line discontinuities in the nematic orientation for C < 0. (a) Paired tubes have 
rims that end on —1/2 defects while (b) closed rims are surrounded by four —1/2 defects (b). 



in a second peak in P(£). This peak position moves to the right as C is increased. We thus observe 
nontrivial changes in the interaction between disclinations due to the geometry of the embedding surface. 
For Cq < 0, there is a proliferation of defects as the inward tubes are formed. The inner side of the 
tubes are decorated with m along the circumference with the tip hosting a pair of +1/2 defects. There 
is a sudden change in the orientation of rh just out side the rim of the tube(see Figj5]). The rim itself is 
thus a line of discontinuity. There are two types of rims, closed ones with four —1/2 defects outside and 
paired open rims ending in two —1/2 defects. 

A variety of thermally excited shapes are also seen. The lowest energy thermal modes among the tubes 
are the spirals, which are always present unless one is at a very low temperature (see figure If). If we raise 
the temperature or decrease the value of ell, the tubes develop branches. More branches are produced 
as the temperature is increased, with every branch having a pair of —1/2 defects at the intersection and a 
pair of +1/2 defects at the tip [23]. Figj6] shows the sequence of snap shots that depict the transformation 
of a sphere with an initial random orientation of the nematic to a tube, when quenched from ell = 
to ell = 3, while keeping k\\ — 5, k± = 0.0, k = 10 and C = 0.4 [30]. The nematic instantly orients on 
the sphere resulting in four +1/2 defects at equal distance from each other. The membrane then begins 
to deform, through formation of protrusions at the defect cores (Figj6ta)), since positive defects favor 
positively curved regions 28 and to satisfy the curvature set by C and K. These membrane protrusions 
are accompanied by production of oppositely charged half defect pairs (Fig |6pD)). The negative defects 
move to the negatively curved intersections and the positive defects move to the positively curved tip of 
the protrusions. When the tip of a protrusion has two +1/2 defects, and its neck two —1/2, that branch 
is stabilized and grows by moving the positive and negative defect pairs away from each other(Figj6jc)). 
Isolated +1/2 and —1/2 defects annihilate(Fig|6jd)). In general two types of branched structures are 
observed, as we increase the directional spontaneous curvature, (a) Broad necked, see FigQg) and (b) 
Narrow necked, see FigJTjh) . 

To analyze the origin of helicity, we neglect the end cap effects of the tube and parame- 
terize it as a canal surface. The total elastic energy, per unit length, is then [24] Wtot = 

~ (zll (l - V^Wf) + ( K+K j~Xr)* ±X2 ) + (iH 1 - C ^ 2x - 1 - HC o r ) + W~ C or)(l - 2x - C£r))Y 
Here r is the radius of the tube and x = cos 2 ip. The energy depends only on the curvature of its 
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FIG. 6: (Color online) The formation of tubes and branches are driven by proliferation of a +1/2 and 

-1/2 defect pairs . See text for details. 



spine curve A, and not on the torsion, hence for A / there is a degeneracy of spiral configurations 
corresponding to different torsions. One can easily see that the energy is at its minimum when A = [24] . 
For a tube of length L, parametrizing the spiral deformations by the extension of the tube along the 
helical axis L z = LsJ\ — Xr, we can estimate the entropy to be of the order fcsTln(^) and for small A 
values the entropy term will dominate the free energy resulting in helical configurations. 

Before we conclude, let us look at how does the model parameters like the spontaneous curvature and 
the density of nematic inclusions compare with experiments. For Cq ~ 1.0, the circumference of the tube 
is about 4 tether lengths. Comparing this with the radius of the lipid tubes obtained in experiments [3T] 
we get the length of the tether to be w 25nm. The value of Cq = 1.0 in real units is thus w (25 nm) _1 , 
which is not far from the suggested value of dynamin intrinsic curvature. Experiments see the dynamin 
rings surrounding the tubes to be made of about 20 units with a pitch of about 15 nm. This will translate 
to about 5 dynamin molecules per vertex in the simulations. 

In conclusion, we have shown that in-plane nematic order couple to curvature, on a deformable surface, 
can lead to non-trivial shapes of vesicles. The deformability of the surface lead to generation of point 
defects and line singularities, which in turn leads to production of tubes and branches. 
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